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INTRODUCTION 


A  number  of  algorithms  have  been  devised  for  the  measurement  of  multifractal 
dimensions.  A  selection  of  such  algorithms  is  described  in  References  1  through  11.  The 
standard  algorithms  yield  the  Hentschell-Procaccia  (ref  1)  fractal  dimension  D^q),  the  f(a)- 
spectrum,  or  other  related  fractal  measures  for  precise  (i.e.,  machine  precision)  data. 

The  consequences  of  noise  on  the  effectiveness  of  fractal  analysis  algorithms  have  been 
addressed  in  References  2  and  3  but  not  in  a  systematic  way.  Reference  2  demonstrated  that  a 
box-counting  fractal  analysis  algorithm  converged  near  "true"  values  for  a  variety  of  Koch 
constructions  and  mappings  based  on  two-  and  three-digit  "pixel-value"  data  for  ^  >  0. 
Reference  3  demonstrated  that  a  generalized  Badii-Politi  algorithm  yielded  the  correlation 
dimension  D{2)  within  about  10  percent  for  a  2*®  point  subset  of  points  on  the  Mackey-Glass 
attractor  having  0.5  percent  random  perturbations. 

This  report  describes  the  results  of  applying  two  multifractal  analysis  techniques  to 
imprecise  fractal  data.  The  correlation  integral  method  (refs  4-6)  and  a  generalized  Badii-Politi 
procedure  are  used  to  determine  the  consequences  of  small  errors  in  the  coordinates  of  points 
comprising  subsets  of  standard  fractal  constructions  (refs  12,13)  in  R^.  The  correlation  integral 
algorithm  is  of  the  customary  form  except  that  automated  procedures  for  selecting  the  scaling 
range  are  employed.  The  generalized  Badii-Politi  algorithm  is  essentially  equivalent  to  that 
applied  by  Kostelich  and  Swinney  (ref  3)  in  their  determinations  of  D{2)',  it  reduces  to  the 
standard  Badii-Politi  algorithm  (ref  7)  as  a  special  case. 


BACKGROUND 

The  Correlation  Integral  Al£orithm 

Analytic  Correlation  Integral  Algorithm 
The  correlation  integral  is  defined  as 

where  Cl  is  the  volume  containing  the  fractal  set,  p{^  is  the  density  of  the  set  at  x,  and  H{x)  is 
the  Heaviside  function. 
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For  a  discrete  fractal  set,  the  correlation  integral,  Eq.  (lb),  takes  the  form 

N  f  .  N  \9-l] 

k  j 


C{q^  =  lim 


1 

q-l 


where  jcj  and  run  over  N-element  fractal  subsets.  It  is  shown  that  the  Hentschell  and 
Procaccia  fractal  dimension  is  then  given  by 

IXq)  = 

£-0*  ln(£) 

The  Numerical  Correlation  Integral  Method 

The  numerical  correlation  integral  method  is  based  on 


1  \  ^ 


\9-l 


£D(q) 


(Ic) 


(2) 


(3) 


where  j  runs  over  the  A^-element  fractal  subset  and  k  runs  over  an  element  randomly  selected 
subset  of  the  fractal  set  in  question.  Equation  (3)  represents  practical  approximations  to  Eqs. 
(Ic)  and  (2). 

The  correlation  integrals  C{q,E)  are  generally  evaluated  for  a  set  of  ^-values  of  interest 
and  a  logarithmically-spaced  set  of  E-values  ehosen  to  cover  the  range  of  scales  in  the  fractal 
subset.  The  smallest  E  is  chosen  to  be  of  the  order  of  the  largest  nearest  neighbor  spacings  in 
the  data  set.  The  largest  E  is  chosen  to  be  approximately  a  diameter  of  the  data  set.  The  fraetal 
dimension  D{q)  is  then  obtained  by  least  squares  fitting  to  the  linear  portion  of  ln[C(^,E)]  versus 
ln[E]. 


An  automated  procedure  was  developed  for  the  selection  of  ranges  of  ln[E]  and  the 
determination  of  D{q)  for  fitting  to  Eq.  (3).  The  limits  on  ln[E]  and  the  values  oi  D{q)  are 
obtained  by  averaging  the  lowest  rms  error  fits  over  ranges  of  consecutive  ln[E].  The  minimum 
number  of  ln[£]  values  is  made  small  enough  to  fit  into  the  linear  range,  yet  large  enough  to 
avoid  small  straight  runs  of  points.  This  procedure  is  described  in  more  detail  in  the  Appendix. 

The  Badii-Politi  Al2orithin 

Analytic  Badii-Politi  Algorithm 

The  Badii-Politi  algorithm  is  designed  for  application  to  discrete  fractal  data.  The 
technique  focuses  on  the  probability  distribution  P{8,n)  of  nearest  neighbor  distances  5  among  n- 
point  subsets.  The  moments  of  the  distribution  of  nearest  neighbor  distances  are  shown  to  vary 
asymptotically  as  where  D{y)  is  a  "dimension,"  which  can  be  simply  related  to  D{q). 
Explicitly, 
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-yWy) 


(4) 


<6^>  H  Af(Y^)  =  f  6'fP(b,n)db  ~  n 
0 

The  dimension  D{y),  which  Reference  3  refers  to  as  "the  dimension  function,"  is  given  by 

DM  =  lim  (5) 

ln(M(Y^)) 

The  Badii-Politi  approach  is  closely  related  to  that  espoused  by  Halsey  et  al.  (ref  14). 

For  example,  Reference  7  demonstrates  how  the  analytic  Ba^i-Politi  algorithm  can  be  applied  to 
determine  the  multifractal  measures  of  Koch  constructions  by  analysis  of  their  generators  along 
the  lines  described  in  Halsey  et  al.  (ref  14).  One  can  translate  the  Badii-Politi  expressions  into 
"standard  notation,"  i.e.,  the  notation  of  Renyi  (ref  15)  or  Halsey  et  al.  (ref  14),  by  the 
replacements: 


y  -  -T  and  D{y)  -  Diq) 

which  imply  that  -y/Z)(y)  ->  q-\  since  t(^)  =  {q-\)D{ci).  Thus,  Eq.  (4)  becomes 

00 

<6-^>  =  M(t,/i)  =  J b-^Pib,n)db  ~  (4’) 

0 

The  Standard  Numerical  Badii-Politi  Algorithm 

Badii  and  Politi  approximate  the  integral  expression  for  M(y,n),  etc.  as 

b]{n)  -  n 


define  the  "y-volume,"  L(y,«)  as 


£(Y,n)  s  n*M(Y,n)  «  ^  6j(n)  ~  =  n’ 

y=i 


where  5  is  taken  as  the  third  nearest  neighbor  distance,  and  compute  D(y)  from  the  scaling  of 
L(y,«)  with  the  number  n  of  randomly  chosen  points  in  the  subsets  in  which  the  third  nearest 
neighbor  distances  are  found.  For  each  n,  a  subset  of  n  points  is  taken  as  a  "reference  set."  (For 
example,  the  average,  defining  L{y,8)  or  M(y,8)  is  based  on  eight  randomly  selected  points.) 
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A  Generalized  Numerical  Badii-Politi  Algorithm 


A  generalized  form  of  the  numerical  Badii-Politi  algorithm  similar  to  that  employed  by 
Kostelich  and  Swinney  (ref  7)  to  determine  the  correlation  dimension  is  used.  It  is  based  on  Eq. 
(4’)  rather  than  Eq.  (5).  The  following  are  the  principal  differences  from  the  numerical 
algorithm  of  Reference  3: 

1.  A  fixed,  relatively  large  reference  set  is  employed  for  all  n.  The  use  of  a  larger 
reference  set,  especially  at  small  n,  substantially  reduces  the  variability  of  M(t^) 
under  different  random  point  selections.  A  substantially  larger  reference  set 
(which  essentially  eliminated  the  dependence  of  results  on  the  random  point 
selection  procedure)  than  that  employed  in  Reference  3  was  employed  in  the 
present  studies. 

2.  The  neighbor  number  77  is  allowed  to  vary.  That  is,  the  scaling  of  nearest 
neighbor  number  distances  for  a  variety  of  neighbor  numbers  (tj)  are  considered. 
Standard  numerical  Badii-Politi  algorithm  results  are  obtained  for  the  special  case 
of  third  nearest  neighbor  (77 =3)  distances. 

Thus,  translating  into  the  notation  of  Halsey  et  al.  (ref  14),  the  generalized  Badii-Politi 
algorithm  takes  the  form 


1 

where  5  can  be  the  distance  to  an  arbitraiy'  77  and  is  independent  of  n. 


(8) 


RESULTS  AND  DISCUSSION 

The  Badii-Politi  and  correlation  integral  numerical  algorithms  were  tested  on  a  variety  of 
constructions  in  R^  for  which  analytic  D{q)  values  are  known.  The  algorithms  were  applied  to 
Euclidean  point  sets,  Koch  asymmetric  snowflake  fractals  based  on  a  four-element  generator 
having  two  0.4  and  two  0.2-length  elements,  symmetric  triadic  snowflakes  (ref  12),  split  snowflake 
halls  (ref  12),  and  the  13-element  generator  Koch  construction  (ref  12).  TTie  attractor  of  the  six¬ 
fold  (D6)  symmetric  chaotic  mapping  described  by  Field  and  Golubitsky  (ref  13)  was  also 
studied.  Imprecise  data  were  produced  by  perturbing  the  x  and  y  coordinates  of  the  fraetal 
points  by  "random"  values. 

The  results  of  extensive  analysis  of  1.97  10“'  point  subsets  of  the  asymmetrie  snowflake 
fractal  are  presented.  Similar  results  were  obtained  for  all  the  fractal  cases  studied.  Reference 
16  demonstrates  that  10’  point  subsets  of  the  asymmetric  snowflake  are  large  enough  to  ensure 
one  percent  convergence  of  the  correlation  integral  method  over  the  range  -25  <  q  <  25. 
Although  we  are  aware  of  no  systematic  study  of  convergence  for  the  Badii-Politi  procedure,  the 
analysis  of  Broggi  et  al.  (ref  17)  suggests  that  1.97  10’  points  are  sufficient  for  Z)(0)  ~  2.  The 
results  of  Kostelich  and  Wnney  (ref  3)  suggest  that  substantially  smaller  point  sets  may  be 
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sufficient.  Also,  Reference  2  established  that  10'^  points  are  sufficient  for  one  percent 
convergence  of  a  box-counting  algorithm  for  the  asymmetric  snowflake  in  the  range  0  <  q  <  25. 

Figure  1  shows  1.23  10^  point  subsets  of  the  asymmetric  snowflake  fractal  with  0.0, 0.5, 

1.0,  and  2.5  percent  random  errors  on  each  axis.  Perturbations  of  the  order  of  0.5  percent  are 
easily  perceived. 

Figure  2  shows  ln[C(-12,£)]  versus  ln[£]  for  1.97  10^  point  subsets  of  the  asymmetric 
snowflake  fractal  with  0.0, 0.05,  and  0.5  percent  random  errors.  The  correlation  integral  values 
were  computed  with  =  N/5  for  40  log-spaced  ^-values  and  were  indistinguishable  from  those 
computed  for  =  N. 

The  vertical  lines  show  the  averaged  upper  and  lower  limits  of  the  range  of  ln[E]  that 
were  selected  by  the  automated  procedure  using  the  twenty  lowest  rms  error  fits  with  runs  of  a 
minimum  of  twenty  consecutive  ln[£]  values.  It  is  apparent  that  they  are  near  "breaks"  in  the 
curves  and  closely  approximate  the  linear  boundaries.  Essentially  the  same  results  were  obtained 
using  the  ten  lowest  rms  error  fits  and  a  minimum  of  ten  consecutive  ln[E]  values. 

These  results  suggest  that  the  0.5  percent  random  errors  introduced  into  the  Mackey- 
Glass  attractor  by  Kostelich  and  Swinney  (ref  3)  (which  are  of  the  order  of  the  uncertainties  in 
their  experimental  Couette-Taylor  data)  were  not  sufficiently  large  to  effect  the  correlation 
integral  analysis  for  the  range  of  hypersphere  radii  (viz.,  -6  <  log2(e)  <  -1)  employed  in  their 
analysis.  One  would  expect  to  find  "breaks"  and  tangible  effects  in  the  correlation  integral 
analysis  for  log2(e)  »  log2(0.005)  «  -7.6,  which  is  substantially  to  the  left  of  the  range  of  log2(e) 
employed  in  their  analyses. 

Note  that  the  linear  portion  for  the  0.5  percent  data  in  Figure  2  comprises  about  one 
order  of  magnitude  variation  in  the  correlation  integral  hypersphere  radii  and  that  this  restricted 
range  of  ln[£']  is  still  sufficiently  large  to  yield  1  percent  agreement  between  the  measured  and 
analytic  D{q)  values  (see  Figure  4b). 

Figure  3  shows  ln[Af(n,y)]  versus  ln[rt]  based  on  the  standard  Badii-Politi  numerical 
algorithm  (tj  =  3)  for  -ylD{y)  +  \  =  TlD{q)  +  \  =  q,  for  q  ~  -12,  ^  ~  0,  and  ^  »  12  for  1.97  10^ 
point  subsets  of  the  asymmetric  snowflake  fractal  with  0.0,  0.05,  and  0.5  percent  random  errors. 
The  Badii-Politi  M{n,y)  values  were  computed  for  twenty  log-spaced  n-values.  The  M(n,y)  values 
shown  were  computed  with  =  NI5.  Unlike  the  data  in  Figure  2  (where  breaks  are  apparent), 
there  is  no  indication  that  the  curves  obtained  from  the  imprecise  fractal  data  should  be  treated 
differently  or  that  analysis  of  the  imprecise  data  will  yield  unreliable  results. 

Figure  4  presents  D{q)  versus  q  curves  determined  by  means  of  the  standard  numerical 
Badii-Politi  algorithm  for  1.97  10^  point  subsets  of  the  asymmetric  snowflake  fractal  with  0.0, 

0.05,  and  0.5  percent  random  errors.  The  solid  line  is  the  analytic  D{q)  versus  q  result.  The 
dense  set  of  points  near  ^  =  0  plus  the  point  near  q  =  5  span  the  range  of  Table  1  in  Reference 
7.  The  standard  Badii-Politi  technique  returns  Z)(-l  <  ^  <  5)  within  1  percent  of  the  analytic 
values  and  D{-5  <  q  <25)  values  within  5  percent  for  precise  fractal  data.  However,  the  results 
are  poor  for  q  outside  this  range  and  small  random  errors  of  the  order  0.05  percent  push  the 
results  up  near  D{q)  =  2.0  (the  random  point  set  value). 
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Figure  5  shows  D{ci)  versus  q  eurves  determined  by  the  eorrelation  integral  algorithm  for 
the  same  data  and  on  the  same  scale  as  Figure  4.  The  solid  line  is  the  analytic  D{q)  versus  q 
result.  The  correlation  integral  results  are  within  1  percent  of  the  analytic  values  of  D{q)  except 
for  a  small  range  near  q  =  -5. 

Figure  6  shows  correlation  integral  algorithm  results  for  D{q)  versus  q  for  1.97  10^  point 
subsets  of  the  asymmetric  snowflake  fractal  with  0.0, 0.05,  and  0.5  percent  random  errors.  The 
solid  line  is  the  analytic  D{q)  versus  q  result.  Although  satisfactory  results  were  obtained  for  1 
percent  fractal  data,  substantial  errors  in  the  measured  D{q)  appear  for  2.5  percent  random 
errors. 


Figure  7  presents  D{q)  versus  q  results  determined  by  the  generalized  Badii-Politi 
numerical  algorithm  using  0.05  percent  data  for  1.97  10^  point  subsets  of  the  asymmetric 
snowflake  fractal.  The  numbers  to  the  left  of  the  curves  designate  the  neighbor  number  used  in 
the  generalized  Badii-Politi  numerical  algorithm.  D{q  >  0)  were  within  2  percent  of  analytic 
values  for  81  <  rj  <  130  and  D{q  <  0)  were  within  5  percent  for  81  <  <  180  for  the 

asymmetric  snowflake.  Similar  effects  were  observed  for  the  other  0.05  percent  constructions 
studied  here. 

Since  the  averaged  81st  nearest  neighbor  distances  are  larger  than  the  0.05  percent 
uncertainty  in  the  fractal  data  set  and  the  averaged  3rd  nearest  neighbor  distances  are  of  the 
order  of  the  0.05  percent  uncertainty  in  the  fractal  data  set,  it  is  not  surprising  that  the  scaling 
for  77  >  81  gives  a  better  representation  than  that  of  the  3rd  neighbors  for  0.05  percent  data. 
(This  effect  of  increasing  77  is  the  motivating  idea  for  the  generalized  Badii-Politi  algorithm.) 

Equivalent  generalized  Badii-Politi  results  were  reported  in  Reference  3  for  the 
correlation  dimension  D(2)  measured  in  Couette-Taylor  flows  for  a  range  of  Reynolds  numbers 
and  in  numerical  realizations  of  the  Mackey-Glass  attractor.  (The  monatonic  decreasing  D(q  > 

0)  with  increasing  77  evident  in  Figure  7  and  for  the  correlation  dimension  in  Reference  3  is  not 
the  general  case  in  the  present  study  or  in  the  results  reported  in  Reference  17.) 

The  \x\[M{n,Y)]  versus  ln[rt]  curves  do  not  reveal  the  range  of  77  which  yield  "best  results." 
For  example,  ln[A/(7i,y)]  versus  ln[«]  curves  that  fit  Eq.  (8)  with  small  rms  errors,  often  yield  poor 
D{q).  Kostelich  and  Swinney  (ref  3)  made  essentially  the  same  comment  with  respect  to  their 
generalized  Badii-Politi  algorithm  determinations  of  the  correlation  dimension. 
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CONCLUSIONS 


The  correlation  integral  method  yielded  better  values  of  D{ci)  than  the  generalized  Badii- 
Politi  technique  in  all  cases  of  noisy  fractal  data  studied  here.  The  numerical  correlation  integral 
method  calculations  reported  here  yielded  D{q)  values  within  1  percent  of  analytic  values  for  0.5 
percent  (and  machine  precision)  representations  of  the  multiffactal  point  sets  studied  except  for 
a  small  range  near  q  =  -5.  {D^q)  values  having  errors  of  the  order  of  one  percent  are  also 
obtained  for  one  percent  representations.) 

The  correlation  integral  method  succeeded  in  extracting  D{q)  from  the  noisy  fractal  data 
because  longer  range  fractal  scaling  is  not  substantially  affected  by  the  imprecision  of  the  data. 
Restricting  the  fitting  to  the  linear  portion  of  the  ln[C(^^)]  versus  ln[£]  curves  excludes  the 
small  "hypersphere"  correlation  integrals  from  the  analysis.  The  range  of  ln[£]  to  exclude  in  the 
determination  of  D{q)  is  apparent  in  the  ln[C(^^)]  versus  ln[£]  plots  and  does  not  require  a 
priori  knowledge  of  the  uncertainties  in  the  data.  The  correlation  integral  method  breaks  down 
when  the  uncertainty  in  the  location  of  the  points  becomes  so  large  that  the  linear  range  in  the 
ln[C(^^)]  versus  ln[£]  plots  is  inadequate. 

The  standard  numerieal  Badii-Politi  algorithm  gave  £>(-1  <  q  <  5)  values  within  1  percent 
of  the  analytie  values  and  D(-5  <  q  <  25)  values  within  5  percent  for  machine  precision  data. 

On  the  other  hand,  the  standard  Badii-Politi  algorithm  failed  at  all  q  when  errors  on  the  order  of 
0.05  percent  were  introduced. 

The  generalized  numerical  Badii-Politi  algorithm  is  less  sensitive  to  noise  than  the 
standard  approaeh.  However,  we  were  unable  to  determine  "optimal"  parameters  by  examining 
ln[A/(T,n)]  versus  ln[n]  curves.  Kostelich  and  Swinney  (ref  3)  and  Broggi  et  al.  (ref  17)  discussed 
the  selection  of  parameters  in  the  Badii-Politi  procedure,  but  were  also  unable  to  establish 
criteria  for  optimal  parameter  selection.  Our  results  and  those  in  Reference  3  imply  that  the 
nearest  neighbor  number  used  in  the  generalized  Badii-Politi  algorithm  must  be  sufficiently  large 
that  the  weighted  average  distances  exceed  the  uncertainty  in  the  positions  of  the  fractal  set. 

D{q  <  0)  for  81  <  Tj  <  180  were  within  5  percent  for  the  asymmetric  snowflake  and 
within  7  percent  for  the  other  constructions.  D{q  ^  0)  for  81  <  rj  ^  130  were  within  2  percent 
for  the  asymmetric  snowflake  and  were  generally  about  5  percent  low  for  the  other  constructions 
studied  here.  This  observation  is  comparable  with  the  results  reported  in  Kostelich  and  Swinney 
(ref  3),  where  300th  neighbor  Badii-Politi  results  for  D(2)  for  the  Mackey-Glass  attractor  is  about 
10  percent  lower  than  the  analytic  value  (i.e.,  6.8). 

We  believe  that  both  the  generalized  Badii-Politi  and  the  correlation  integral  algorithm 
can  be  successfully  applied  to  impreeise  data  in  high  D  eases.  However,  problems  such  as 
restricted  linear  scaling  ranges  and  lack  of  convergence  remain  unresolved. 

A  final  point  in  regard  to  analysis  of  impreeise  fractal  data  is  in  order.  If  the  fractal 
point  set  in  question  is  the  result  of  invariant  embedding  of  a  chaotie  time  series,  then 
preprocessing  of  the  data  along  the  lines  described  by  Kostelich  et  al.  (ref  18)  may  affect 
significant  noise  reduction. 


7 


REFERENCES 


1.  H.G.E.  Hentschel  and  I.  Procaccia,  Physica  D,  Vol.  8,  1983,  p.  435. 

2.  L.V.  Meisel,  M.  Johnson,  and  P.J.  Cote,  Phys.  Rev.  A,  Vol.  45,  1992,  p.  6989. 

3.  E.J.  Kostelich  and  H.L.  Swinney,  Physica  Scripta,  Vol.  40,  1989,  p.  436. 

4.  K.  Pawelzik  and  H.G.  Schuster,  Phys.  Rev.  A,  Vol.  35,  1987,  p.  481. 

5.  H.  Atnianspacher,  H.  Scheingraber,  and  G.  Wiedenmann,  Phys.  Rev.  A,  Vol.  40,  1989,  p. 
3954. 

6.  G.  Paladin  and  A.  Vulpiano,  Lett.  Nuovo  Cimento,  Vol.  41,  1984,  p.  82. 

7.  Remo  Badii  and  Antonio  Politi,  J.  Stat.  Phys.,  Vol.  40,  1985,  p.  725. 

8.  A.  Block,  W.  von  Bloh,  and  H.J.  Schellnhuber,  Phys.  Rev.  A,  Vol.  42,  1990,  p.  1869. 

9.  X.J.  Hou,  R.G.  Gilmore,  G.B.  Mindlin,  and  H.G.  Solari,  Phys.  Lett.  A,  Vol.  151,  1990,  p. 
43. 

10.  Mogens  H.  Jensen,  Leo  P.  Kadanoff,  Albert  Libchaber,  Itamar  Procaccia,  and  Joel 
Stavans,  Phys.  Rev.  Lett.,  Vol.  55,  1985,  p.  2798. 

11.  L.S.  Liebovitch  and  T.  Toth,  Phys.  Lett.  A,  Vol.  141,  1989,  p.  386. 

12.  For  example,  B.B.  Mandelbrot,  Fractal  Geometry  of  Nature,  Freeman,  New  York,  p.  1983. 
(13  element  construction:  p.  69,  split  snowflake  halls:  p.  146). 

13.  M.  Field  and  M.  Golubitsky,  Computers  in  Physics,  Vol.  4,  1990,  p.  470. 

14.  Thomas  C.  Halsey,  Mogens  H.  Jensen,  Leo  P.  Kadanoff,  Itamar  Procaccia,  and  Boris 
I.  Shraiman,  Phys.  Rev.  A,  Vol.  33,  1986,  p.  1141. 

15.  A.  Renyi,  Probability  Theory,  North  Holland,  Amsterdam,  1970. 

16.  L.V.  Meisel  and  M.A.  Johnson,  "Convergence  of  Numerical  Box-Counting  and 
Correlation  Integral  Multifractal  Analysis  Techniques,"  submitted  to  Computers  in  Physics. 

17.  G.  Broggi,  J.  Opt.  Soc.  Am  B,  Vol.  5,  1988,  p.  1020.  M.  Ravani,  B.  Derighetti,  G.  Broggj, 
and  E.  Brun,  J.  Opt.  Soc.  Am  B,  Vol.  5,  1988,  p.  1029. 

18.  E.J.  Kostelich  and  J.A.  Yorke,  Physica  D,  Vol.  41,  1990,  183.  E.J.  Kostelich,  Physica  D, 
Vol.  58,  1992,  p.  138.  E.J.  Kostelich  and  T.  Schreiber,  Phys.  Rev.  E,  Vol.  48,  1993,  p. 

1752. 


8 


Figure  1.  1.23  10'*  point  subsets  of  the  asymmetric  snowflake  fractal  with 
(a)  0.0,  (b)  0.5,  (c)  1.0,  and  (d)  2.5  percent  random  errors. 


ln[E] 


Figure  2.  ln[C(-12^)]  versus  ln[iE]  for  1.97  10^  point  subsets  of  the  asymmetric 
snowflake  fractal  with  0.0  (solid  line),  0.05  (filled  circles),  and  0.5 
percent  (x’s)  random  errors.  The  vertical  lines  are  the  averaged 
upper  and  averaged  lower  limits  automatically  selected  for  the  range 
of  ln[E].  The  limits  for  the  precise  data  are  unmarked,  for  0.05 
percent  data  are  marked  with  a  filled  circle,  and  those  for  0.5 
percent  data  are  marked  with  an  x. 
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O  □  □ 


□  □  □ 


Figure  6.  D{q)  versus  q  results  using  the  numerical  correlation  integral  algorithm 
for  1.97  10“'  point  subsets  of  the  asymmetric  snowflake  fractal  with  0.0 
(open  circles),  0.5  (filled  circles),  1.0  (x’s)  and  2.5  percent  (open 
square)  random  errors.  The  solid  line  is  the  analytic  solution. 


APPENDIX 


Automated  procedure  for  correlation  integral  computation  of  D{q): 

1.  Specify  two  numbers: 

a.  A  minimum  number  of  consecutive  logarithmically-spaced  ln[£]-values  to  be 
contained  in  an  allowable  "linear  range."  This  number  should  be  large  enough  to 
avoid  local  straight  runs  of  points  and  small  enough  to  fit  into  the  "linear  range." 

b.  The  number  of  ranges  to  include  in  determining  D{q).  Since  the  \n[C{q^)]  versus 
ln[£]  points  tend  to  oscillate  around  a  straight  line,  "best"  results  are  obtained  by 
averaging  over  sets  of  consecutive  ln[£]-values  that  terminate  at  high  and  low 
ln[C(^,£)]  values. 

2.  Make  least  squares  fits  to  Eq.  (3)  for  all  ranges  of  ln[£]  consistent  with  l.a. 

3.  Sort  the  results  on  rms  errors  in  slope. 

4.  Average  values  over  the  lowest  rms  error  cases  up  to  the  number  specified  in  l.b.  and 
return  the  averaged  smallest  ln[E]  value,  the  averaged  largest  ln[£']  value,  and  the 
averaged  D{q)  value. 
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